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Abstract 



Oh' We calculate the mass splittings within the heavy baryon isospin multiplets X c (b) and S' c ^ in 

yi^ . chiral perturbation theory to leading one-loop order. The pattern of the mass splittings in 

' the E c iso-triplet, which is different from that of any other known isospin multiplet, can be 

,JZ< , explained. We predict m~i+ — rano = —0.2 ± 0.6 MeV, m s m — m-i- = —4.0 ± 1.9 MeV and 

the mass of the to be 5810.3 ±1.9 MeV. 
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1 Introduction 

> 

•^H , 

Mass splittings within isospin multiplets of hadrons appear due to both the mass difference between 
the u and d quarks and electromagnetic (em) effects. Since the d quark is heavier than the u, usually 
the hadron with more d quarks is heavier within one isospin multiplet. For instance, the neutron 
(udd) is heavier than the proton (uud), and the K°(ds) is heavier than the K + (us). There is only 
one exception to this pattern, the S c iso-triplet, consisting of the £+ + (cu-u), the £+(cu<i) and the 
YP c (cdd). The mass splittings within the S c iso-triplet are measured [I] 

Ai c = - m s o = -0.9 ± 0.4 MeV, 

A 2c = m E ++ - m E o = 0.27 ±0.11 MeV. (1) 

Remarkably, the state with two u quarks has the largest and the one with a u and a d quark has 
the smallest mass. 

Only recently some of the bottom cousins of the S c , the E^, were observed by the CDF Col- 
laboration [2]. Their masses are for the buu state m s + = 5807.8 ± 2.7 MeV and for the bdd state 
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m s - = 5815.2 db 2.0 MeV, respectively — their neutral partner Ejj has not been observed yet. 

b 

Thus, here the natural ordering of the states seems to be restored. On the other hand, heavy 
quark symmetry relates baryons containing a b quark to those with a c quark, which makes this 
different pattern even more puzzling. In this work we investigate the origin of these patterns, to- 
gether with those in the S' c / fc s doublets, using chiral perturbation theory (CHPT) and heavy quark 
symmetry. In this way we can include both sources of isospin violation in a way consistent with 
QCD |3j. Our work is a straightforward extension of analogous studies for the nucleon [U [5] (for 
pioneering studies, see [3[7])- 

The isospin splittings for heavy baryons were already studied in various quark models in Refs. [HJ 
El [TUl IH 121 021 [HI [151 HE] • The advantages of our investigation are that (i) , since we use an 
effective field theory, the theoretical uncertainty can be estimated, (ii) for the first time meson loop 
corrections are considered — they turn out to be numerically significant for the charm baryons, 
and (Hi) this study investigates all isospin splittings at the same time. We compare our results to 
those of the quark models below. 

Our paper is organized as follows. In Section [2l we construct the SU(3) chiral effective La- 
grangians responsible for the mass corrections of the sextet heavy baryons at order 0{p 2 ). Here, 
p denotes the expansion parameter of the underlying effective field theory. The strong Lagrangian 
is proportional to the quark mass and the em Lagrangian is constructed including virtual photons. 
As we will show, the operator structure of the effective field theory for the heavy-light baryons is 
richer than the one of the light quark sector since heavy and light quarks must be treated differ- 
ently. Calculations up to 0(p 3 ) are performed in Section [3l A brief summary is given in the last 
section. Some technicalities are relegated to the appendices. 

2 The chiral effective Lagrangians 

In CHPT the quark mass difference enters explicitly through the quark mass matrix. The inclusion 
of the virtual photons has been first considered systematically for the three-flavor case in Ref . [T7j . 
Chiral Lagrangians with virtual photons have been constructed for the study of isospin symmetry 
breaking phenomena in mesons and baryons with light up and down quarks (see, e.g., [5l [T8 l [T9| 
GOlEIlEaEaiMlESliaEB^ incomplete list). Recently, this technique was used to 

study the interaction between Goldstone bosons and heavy-light mesons and the isospin breaking 
decay width of the L>* (2317) [30]. 

In this paper, we use the technique of SU(2) chiral perturbation theory to study the mass 
splittings within the heavy baryon isospin multiplets. The u, d quark masses and the electric 
charge e are counted as small quantities. They are booked as m u ,irid ~ 0(p 2 ), and e ~ 0(p) 
as usual, where p denotes a small momentum with respect to the typical hadronic scale of about 
1 GeV. Both types of isospin-violating effects ~ (m u — m^) and ~ e are taken into account here 
in a systematic manner to the order 0(p 3 ). These effects can be accounted for in three-flavor 
CHPT to study the mass splittings within the SU(3) multiplets systematically. However, similar 
to the case of CHPT for light baryons, the SU(3) breaking contributions from the kaon-baryon 
and eta-baryon loops are large, see e.g. [31], [321 [33], which makes the convergence of the chiral 
expansion problematic. We will therefore treat each heavy baryon isospin multiplet separately to 
0(p 3 ) in two-flavor CHPT, and relate the low-energy constants (LECs) to 0(p) and 0{p 2 ) through 
SU(3) relations (see also the discussion in Section [3]). On the other hand, the splittings within the 
multiplets are well behaved. This procedure is equivalent to starting from the SU(3) Lagrangian 
and to calculate only the leading order SU(2) loops, i.e., the pion-baryon loops. 

In this section, we will construct the SU(3) chiral effective Lagrangians pertinent to the mass 
corrections of the heavy baryons to order 0(p 3 ) (for similar works considering the strange quark 
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as heavy, see Refs. [SI ESI ESI EI])- 

In order to construct the chiral effective Lagrangians, the following building blocks are necessary 
(we employ the standard nonlinear realization of chiral symmetry) 

TT (V2iA 2 u 

U = exp I I , u A = U, 



in 



VJJu\ 



u 
1 



X+ = u ] xu ] + uxu, 



Q± = | [^Qu±uQ^) , (2) 



where F w is the pion decay constant^, <fi collects the Goldstone boson fields 



\ 




~^° + j E V K» |. (,) 
K- K° 

and is the photon (em) field. The diagonal quark mass matrix and the charge matrix are 

X = 2B ■ diag {m„, m d , m s } , 

Q = e- diag {2/3, -1/3, -1/3} , (4) 

in terms of Bq = \(0\qq\Q)\/F% and the elementary electric charge e (e > 0). Under SU(3)lxSU(3)#, 
Uyrj, x+ an d Q± transform as 

O -► hOh\ (5) 

where the compensator field h is an element of the conserved vector subgroup SU(3)y. 

The E c iso-triplet and the Ef c iso-doublet belong to the symmetric sextet and the A+ and the H c 
doublet belong to the anti-symmetric triplet in the flavor SU(3) classification. We use the following 
matrix representation in accordance with the notation of Refs. |38|. [39] 

V2S++ £+ ~'+ \ / A+ 

V2^ c E' c ° , B- 3c = -A+ E° c | . (6) 

k° >mJ v-- c + -sg o 




Under SU(3)l xSU(3)i?, the transformation laws of the charmed baryon fields are 

5 6c -> ZlSg^ 7 , fig, -> fcBg c /l T , (7) 

where /i T is the transpose of /i. The matrices for the bottom baryons can be obtained replacing c 
by b and decreasing the electric charge of every state by one unit. 

Let -BgQ, Bqq (Q = c,b) denote the heavy baryon fields, and m 3 Q, rh§Q their masses in the 
chiral limit, respectively. Analogous to the effective chiral Lagrangian for the pion-nucleon system, 
the lowest order Lagrangian involving the sextet and anti-triplet heavy baryon fields can be written 
as [381 HQ] 



= \ (B- 3Q (i^D» - m- 3Q )B- 3Q ) + (Bsq^D* - m 6Q )B 6Q 



2 

1 - 1 - 1 

+ {BqqW^Bqq) + ~g 2 (BqqIhI^Bzq) + h.c. + -g 3 {B^^u^B^q) . (8) 



1 Strictly speaking, this should be the pion decay constant in the chiral limit. To the accuracy we are working, 
however, we do not need to differentiate this from its physical value. 
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with 

Dfi = Dfj,- iQb+A^, u^ = u^- 2Q B -A^ . (9) 
The chiral covariant derivative on the baryon fields is given by 

r M = i (vta^u + ud^ . (io) 

The charge matrix of the heavy baryons Q b , which gives the correct minimal coupling of the heavy 
baryons to photons, is constructed as 

„ „ I e • diag {2, 0, 0) , for the charm baryons, . , 

Q B = 2Q + q h I = { r ( U ) 

I e • diagjl, —1, —1} , for the bottom baryons, 

where q^ is the charge of the heavy quark, and I is a 3 x 3 unit matrix. Qb± is defined as 

Qb± = ^ (u^Q B u ± uQbu^ . 

At 0(p 2 ) the strong Lagrangian pertinent to the corrections of the masses is given by the terms 
containing one power of the quark mass matrix 

42. = - (Bq (aix+ + « 2 (x+» ^Q> • (12) 

The em Lagrangian at 0(p 2 ), parameterizing hard virtual photons, is more complicated than the 
one of the pion-nucleon system [HE]. The terms, which contribute to the mass corrections of the 
sextet heavy baryons and are quadratic in the light quark charge matrix read 

6$ Q = -Fl (B 6Q [A, {Ql - Ql) + t3 x Q + (Q+) + 2 (Q\ - Ql) + 03 (Q 2 + + Q 2 >] B 6Q ) 

-F|/3 4 (QIb 6 qQ + B 6 q) . (13) 

The em Lagrangian given above only deals with the hard photons exchanged between the light 
quarks in the heavy baryons. In addition we need to add terms that parameterize the em inter- 
actions between the heavy and a light quark. Since the heavy quark can be viewed as static, its 
charge acts as a static background field that transforms as a scalar under SU(3)lxSU(3)r. It 
is easy to see that in this way we get one additional independent structure that contributes to the 
isospin splittings, namely 

4?1 = 4o - F ^ lh (B6 Q Q + (q h l)B 6Q ) . (14) 

As we will show below, it is this term which makes the mass splitting pattern within the S c iso- 
triplet different from that within the E& iso-triplet. The analogous mechanism is also present in 
the mentioned quark model calculations. 



3 Mass splittings within the heavy baryon isospin multiplets 

As mentioned above, our strategy is to relate the LECs for different heavy baryon isospin multiplets 
using SU(3) relations, but, since we are only after isospin splittings, calculate only the pion-baryon 
loop corrections. If we start from the SU(2) Lagrangians for each of the isospin multiplets, the 
chiral limit masses of the Sq and E'q are different, and the difference contributes to the SU(3) 
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mass splittings of 0(p 3 ). It is important to stress that the differences of the LECs between their 
SU(2) and SU(3) values starts contributing at order 0(p 4 ) and thus is beyond the accuracy of our 
calculation!! Therefore, what we do is equivalent to use the SU(3) Lagrangians given in Section [2] 
to do SU(2) calculations. Indeed, the number of independent operators matches for the two cases, 
since the contribution of the (3\ term to the baryon masses vanishes for SU(3) because the three- 
flavor quark charge matrix is traceless, and the (5q term for SU(3) can be translated into the (3\ 
term for SU(2) by using the Cayley-Hamilton relation for 2x2 matrices 

Ql-Q- = Q+ (Q+) + \ ((Ql - Q-) - (Q+) 2 ) • (15) 
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3.1 Mass splittings to 0(p 2 



At leading order, there is no mass splitting within isospin multiplets. However, at next-to-leading 
order the terms in the 0(p 2 ) Lagrangians generate mass splittings. From Eqs. (|12ll3p . we get 

A i? = (™s+ - ™Eg) = 2 ai B (m u - m d ) + -F 2 e 2 ((3 Q - 2(3 4 + 6(3 lh ), 
A< 2c = ( m s++ ~ m s°) ( ) = 4ai-B (m n - m d ) + ^F 2 e 2 (f3 + /3 4 + 6f3 lh ), 
A 3? = ( m ~>+ ~ ™h?) = 2 ai B Q (m u - m d ) + -F 2 e 2 (f3 - 2/? 4 + 6fr h ) = Ag>, (16) 
for the charm baryons S c and r! c . Similarly, for the bottom baryons we have 

A ?b = ~ = 2ai J B (m u - m d ) + -F 2 e 2 ((3 - 2/? 4 - 3ft h ), 



A 26 = ( m s+ ~ m s- ) = ^iB (m u - m d ) + -F n e (f3 + /? 4 - 3/%), 

(^) 1 

A 3? = ( m s' b o ~ "V r ) = 2ai J B (m u - m d ) + -F 2 e 2 (p ~ 2p4 - Wlh) = ■ (17) 

Note that the a\ and (3q terms always appear in the same linear combination, which means that 
the strong contribution, the a\ term, cannot be disentangled from the em contribution without 
additional information. This is completely analogous to the case of the neutron-proton mass 
splitting, see e.g. Ref. [3]. The (5\h term has a different sign for charm baryons and bottom baryons, 
hence it is expected to induce a different interference pattern of the various contributions. 
From the relations given above we find 



m s t+ — m^,oj — — m £0 j = 0(p 3 )> 

(m s ,o - rn B ,-j - (m s o - m E -J = 0(p 3 ), 

+ m s - — 2ra s oj — ^m s ++ + m s o — 2m E +^ = 0(p 3 ). (18) 
The last relation is obtained invoking heavy quark symmetry. 

3.2 Mass splittings to 0(p 3 ) 

The first non- vanishing loop corrections to the baryon masses appear at order C(p 3 ). At this 
order formally both photon loops as well as pion-baryon loops contribute. Remarkably, QCD does 



2 The precise matching between two- and three-flavor versions of CHPT is discussed in Refs. |41l 1421 [43] . 
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Baryons 


Loops 








S++7T 




h , A+7T+ 




S++7T- 






£° 


S+7T-, 


sv, 


A+vr- 




S'+TT 


S / C 7T+ 


E+TT° E°TT+ 




E'+TT" 




~°7T 



Table 1: Pion-baryon loops contributing to the charm baryon mass corrections. 



not allow for a counterterm at this order and consequently the 0(p 5 ) pieces of the these loops 
are finite. We start with the latter kind of loops that are to be constructed from two vertices of 
0(p). A complete list of loops is given in Table [JJ for the charm baryons considered here. Since 

71 7l + 71 + 

S N S ^ / ^ 

/V /V /V 

/ \ / \ / \ 
I 1 I I I I 

S ++ E ++ E ++ E ++ E * E ++ E ++ A + E ++ 

c c c c c c c c c 

(a) (b) (c) 

IX ' 71 + 71° 

^ ^ ^ - - 

/ N / V / N 

/ \ / v / \ 
I 1 I 1 I 1 

E + E ++ E + E + E E + E + A + E + 

C CCC CcC cc 

(d) (e) (f) 

71 71' 71' 

.'«-s 

/ V / V / N 
I \ I \ / \ 
I I I I I I 



E E E E E + E E A + E 

C CCC CCC c c 

(g) (h) (i) 

Figure 1: The pion-baryon loops contributing to the self-energies of the E+ + (a, b, c), E+ (d, e, 
f) and S° (g, h, i). The black dots in the charged pion propagators denote the electromagnetic 
insertions at 0(p 2 ). 



in the power counting the pion mass difference is of the same order as the pion mass itself, in 
the loops we are to use the physical pion masses (for a detailed discussion of this point, see e.g. 
Ref. |26j). On the other hand, to the order we are working, the masses to be used for the baryons 
of the same isospin multiplet are the same. Therefore, the leading SU(2) loop contributions to 
the mass corrections of the E'^~ and E'® as well as S+ + and £j? are equal, as can be seen from 
Table [TJ Hence there is no loop correction for the corresponding mass differences at 0(p 3 ). The 
pion-baryon loops for the S c self-energies are shown in Fig. [TJ Contrary to the case of the nucleon 
mass differences, here it is not straightforward to use the heavy baryon formalism to calculate the 
pion-baryon loops, since the pion-A c contribution generates a cut. It is thus more convenient for 
us to evaluate the integrals using the covariant method of infrared regularization as derived by 
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Figure 2: The photon-baryon loop. 



Becher and Leutwyler [33]. Some remarks on the method and the relevant integrals are given in 
Appendix [Al 

The photon-baryon loops are shown in Fig. [2j Formally they also contribute at 0(p 3 ). However, 
it can be shown that they vanish to this order. Since the baryon mass in the loop is equal to the 
mass of the external legs, for this loop we may use the integral representation of the heavy baryon 
formalism. Then the vanishing of the loop follows from the absence of a mass scale in the integral. 
This result also holds, when the infrared regularization is employed, as outlined in Appendix [Al 

Therefore the mass splittings for the charm baryons E c and r! c are to third order in the chiral 
expansion 

m s + - m 2 o = Ag> + A[° c op (mx c ,m Ac ) + 0(p 4 ) , (19) 
m Et+ - m s o = Ag + 0( P 4 ) , (20) 
m s ,+ - m B ,o = Af] + 0(p A ) . (21) 
The mass splittings for the bottom baryons and r! b are 

(22) 
(23) 
(24) 

where the results to 0(p 2 ), the A^ 2 ) were given in the previous subsection. As derived in Ap- 
pendix [Bl the loop functions to 0(p s ) are given by 

^K.^) = (Ml -Ml) 

+ReS( c )(m SQ ;M ,m AQ ) -ReS( c )(m SQ ;M ± ,m AQ ) , (25) 

with Mq (M±) the neutral (charged) pion mass. The explicit expressions for the loop functions are 
given in Appendix [Bl 

3.3 Numerical results 

When calculating the loops, we take the physical values for the masses 

M± = 139.57 MeV, M = 134.98 MeV, 
m Sc = 2453.56 MeV, m Ac = 2286.46 MeV, 

m S(j = 5811.5 MeV, m Ab = 5620.2 MeV. (26) 
The values of g\ and gi can be estimated based on SU(6) [38] 

gi = \gf, to = -\J\gf, (27) 





— m v - 

^b 


= A (2) 


+ Ai° op (m Eb ,m A J + 0(p 4 




— m v - 
^b 




+ <V) , 


mnfl 


— m~i- 


= A (2) 


+ o( P 4 ) , 
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where is the coupling constant for the single quark transition u —>■ d. A value of g u J^ = 0.76 
gives the correct nucleon axial coupling constant gA = (5/3)^ rf = 1.27, correspondingly g\ = 1.02. 
It also gives g^ = —0.62, the absolute value of which is close to the empirical one \g^ = 0.58 ±0.04 
obtained from the measured decay width r(£+ + —* A+7r + ) = 2.23 ±0.30 MeV. The only unknown 
parameters are the LECs in the 0{p 2 ) La grangians. There are effectively three 

7 = 2aiB (m u - m d ) + ^e 2 F 2 (3 , 
fa = e 2 F 2 fa, 

fa h = e 2 F 2 (3 lh . (28) 

Totally there are four known isospin mass splittings of the sextet heavy baryonsjl 

m s + - m s o = -0.9 ± 0.4 MeV, 
m s ++ - m s o = 0.27 ± 0.11 MeV, 
m + - m = -7.4 ± 3.4 MeV, 

b b 

m H ,+ - mg/o = -2.3 ± 4.2 MeV. (29) 

The first three will be taken to determine the LECs 7, (3^ and (3\h because they have the smallest 
uncertainty. 

Using the physical value of the pion decay constant F w = 92.4 MeV, we get the contribution of 
the loops to m s + — m s o 

A l °° p = (-0.32 ± 0.15) + (-0.41 ± 0.06) MeV 

= -0.73 ±0.16 MeV, (30) 

where in the first line the numbers in the first parenthesis are from the tt-T, c loops, and those in the 
second parenthesis from the n-A c loops. The uncertainty of the loops is in general controlled by 
the expansion parameter of CHPT, namely x = M/h x , with the chiral symmetry breaking scale 
A x ~ 1 GeV. Since we do not really know the value of g%, which enters in the tt-T, c loops, for those 
we estimate the uncertainty conservatively as being of the order 4%. Since the coupling constant of 
the S c to A c and tt can be extracted from experiment, for the second contribution we use directly 
the uncertainty that results from that extraction (see above). Note that the given uncertainty is 
at the same time numerically of order x- Thus the uncertainty estimate is consistent with what is 
expected from the chiral expansion. The LECs are then determined as 

7 = -2.5 ± 1.1 MeV, 
(3 4 = 0.6 ± 0.9 MeV, 
(3i h = 2.6 ± 1.1 MeV. (31) 



The two different contributions in Eq. (|30p are comparable, and they are considerably smaller 
than the individual strong and em contributions at 0(p 2 ), see Eqs. (|16p and (I3ip . showing good 
convergence of the chiral expansion. The mass splitting within the S' c doublet can be predicted 

m H ,+ - m E ,o = m s + - m E o - A^" p = -0.2 ± 0.4(exp) ± 0.4(th) MeV, (32) 

where the first uncertainty coming from the uncertainty of m-^+ — m E o is experimental, and the 
second one is theoretical. It comes from neglecting the 0(p 4 ) contribution and was estimated 
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£ c ++ - £° 


l— 'c 


£° - 


K - ^ 


"6° a 6 


Exp. n 


-0.9 ±0.4 


0.27 ±0.11 


-2.3 ±4.2 




-7.4 ±3.4 




Our work 


-0.9 ±0.4* 


0.27 ±0.11* 


-0.2 ±0.6 


-4.9 ± 1.9 


-7.4 ± 3.4* 


-4.0 ± 1.9 


m 


-0.9 ±0.4* 


0.27 ±0.11* 




-6.9 ± 1.1 


-7.4 ±2.3* 




m 


-0.83 


-0.01 


-1.48 








mi 


-0.73 


0.28 


-3.20 


-3.95 


-6.12 


-6.16 


ma 


-0.2 


1.4 




-3.7 


-5.6 




rrm 

P2] 


-0.5 


3.0 


-1.0 


-5.6 


-7.1 




53] 


-0.7 


0.5 


-1.2 










-0.40 


0.84 










na 


-0.36 


1.20 


-0.30 


-2.51 


-3.57 




usi 


-0.33 


0.37 


-0.20 









Table 2: Comparison of our results with the experimental data and the results from quark models 
(units are MeV). The numbers marked by * are used as inputs. 



by taking one half of the leading loop contribution. Because the width of the £& has not been 
measured so far, g<i for the bottom baryons cannot be determined from the data. Taking the same 
values as before for g\ and 52 > as dictated by heavy quark symmetry, we get the mass difference 
between E° and S^" at C(p 3 ) 

A l °° p = (-0.3 ± 0.2) ± (-0.6 ± 0.1) MeV 

= -0.9 ± 0.2 MeV. (33) 

where in the first line the numbers in the first and the second parentheses are from the n-'Eb 
loops and the n-Ab loops, respectively. The uncertainties were estimated as in case of the charm 
baryons. According to Eq. (|17p . the mass of the SP, which has not been measured yet, is predicted 
to be 

"iso = 5 ( m s+ + m s- - + A 1 ^ 015 

= 5810.3 ± 1.8(exp) ± 0.5(th) MeV. (34) 

We can also predict the mass difference between the E' b ° and E' b ~ 

m = /o — to-/- = — (m^+ — m v - — 3a) 

-fc 2 V H J 

= -4.0 ± 1.8(exp) ± 0.5(th) MeV. (35) 

In Table [2j the results of our work are summarized and a comparison with those obtained in quark 
models is given. 

3 The first two mass splittings are given in PDG pQ. The last two mass splittings are evaluated by taking the 
difference of the masses given in PDG: m E+ = 5807.8±2.7 MeV, m £ - = 5815.2±2.0 MeV, m s ,+ = 2575.7±3.1 MeV, 
and m s ,o = 2578.0 ± 2.9 MeV. 
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4 Discussions and summary 



In this work, we have calculated the mass splittings within the heavy baryon isospin multiplets 
£ c (&) and to 0(p s ) in the chiral expansion. Our main results are given in Eqs. (|19H24j) and 
in Table [21 To arrive at these results, we constructed both the strong and the em Lagrangians 
at 0(p 2 ) which are responsible for the mass corrections. In contrast to mass splittings in light 
quark baryon multiplets, there is an additional operator that describes the hard virtual photons 
exchanged between the heavy quark and light quarks accompanied by a LEC (3\h- Remarkably, 
this term has a different sign for the charm baryons and the bottom baryons. This is due to the 
fact that the sign of the electric charge of the charm quark is different from that of the bottom 
quark. It is the different interference between this term and the other terms that drives the 
mass splittings within the E c iso-triplet to have a different pattern compared to any other known 
isospin multiplet. This leads one to expect that the isospin mass splittings in the charm hadrons 
are always different from those in the bottom hadrons even if the heavy quark symmetry were 
exact. Besides the heavy baryons considered in this paper, the D and 5-meson mass splittings, 
m D ± — m D o = 4.78 ± 0.10 MeV and m B o — m B ± = 0.37 db 0.24 MeV p] are a nice example for the 
effect, although the ordering does not get changed here. 

There is no loop contribution to the mass splitting between the two Ef c baryons, and we predict 
m~/+ — m^io = —0.2 ±0.6 MeV. The present data for the masses of the S' baryons are not accurate 
enough yet to test this prediction. For the S& states, the j3\h term interferes constructively with 
the other terms and hence the loop corrections are less important. The mass of the E? and the 
mass difference m = /o —m~i- are predicted to be 5810. 3±1. 9 MeV and — 4.0±1.9 MeV, respectively, 
which can be tested in future experiments. 
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A Infrared regular izat ion 

The finite masses of the baryons in the chiral limit spoil the correspondence between the loop 
expansion and chiral expansion if one uses conventional dimensional regularization [45]. Here we 
follow the infrared regularization (IR) method developed by Becher and Leutwyler [44J to overcome 
this problem. The IR method has been extended to two loops [46] (see also Ref. |47|). and to the 
cases with spin-3/2 fields [IB] and spin-1 fields [Ml [50] (for reviews, see [51| [52] ) . 
In the IR method, the scalar loop integral 

I r d d k 1 

~ % J (2n) d (M 2 - k 2 ) [m 2 - (P - k) 2 ] ' ( 

with M and m being the masses of the pion and baryon, respectively, is split into two parts, one 
being infrared singular and the other being infrared regular, H = I + R. Only the singular part /, 
which is of order 0(p), makes the expansion in loops to coincide with the chiral expansion, hence 
leads to a consistent power counting. The regular part R can be expanded in polynomials in M, 
hence it can be absorbed into the LECs order by order. 
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In any regularization such as the IR which has a consistent power counting for loops, the 
photon-baryon loops with each of the vertices being of the 0(e) order should be counted as 



Q(e p) = 0{p ). In the photon-baryon loops, when taking P 



m 



which is necessary for 



calculating chiral corrections to the baryon mass, there is no quantity of the order of 0(p) since 
the photon is massless. Therefore, the photon-baryon loops should vanish for calculating the mass 
shifts of baryons. Such an argument is supported by explicit calculations in the IR method. 
The self-energy from any of the pion-baryon diagrams has the form 



-i^ n (P;M,m) 



d d k 



i(P- ft + m) 



AF 2 J (2vr) d ' 10 (k 2 — M 2 + ie) \(P - k) 2 - m 2 + ie] 



(A.2) 



where (n) is a diagram label, P is the external momentum, and (i = 1,2) are the £q£q7t 
and SqAq7t (Q = c,b) coupling constants of the lowest order Lagrangian, Eq. (|8|). After a few 
manipulations one gets 



T^ Q {P-M,m) = J>j-(P + m) M 2 I(P 2 ) + (m- P) PI^(P 2 ) - A Sf 
In the IR method, the single baryon loop 

d d k 1 



(A.3) 



(2vr) d k 2 - m 2 



vanishes. The expressions for the loop functions I(P 2 ) and I^\P 2 ) are given in Appendix [Bl 

Now let us focus on the Taking P = m s ++, the expressions for the diagrams Fig. Q]Ja 

and (b) are very simple because the term proportional to (m— P) does not contribute, 

£( a W;M ,m E J- ^ 



32irF 2 ' 



(A.4) 



where the chiral limit mass mj; e in the loop has been replaced by the physical mass since the 
contribution of the difference is of higher order. Note that up to 0(p 3 ), we do not need to 
distinguish the masses of baryons with different electric charges in loops, so that ms c , instead of 
m s ++, is used for the arguments of the loop function. Replacing the neutral pion mass Mq by 
the charged pion mass M±, the expression for diagram (b) is obtained. The contributions of the 
diagrams (a) and (b) do not depend on the baryon mass. 

The expression for diagram (c) is much more complicated since the term proportional to (m— P) 
has a finite contribution 



£( c )(m Ec ;M ± , 



A 



(m Sc + m Ac ) M|l(m! ) + (m Ac - m s Jm Sc I (1) (ml) 



(A.5) 



where the loop functions I(P 2 ),I^(P 2 ) are replaced by their finite parts I(P 2 ), 1^ (P 2 ) (sub- 
tracting the A parts), and the divergences can be absorbed in the counterterms at 0(p 4: ) which 
are not considered here. The expressions of the infrared singular parts of the loop integrals are 
given in Appendix [Bj Here we use the expansion of the I(P 2 ) up to 0(p). The chiral expansion 
of I{P 2 ) up to 0(p), accounting for the cut due to the opening of the A c 7r channel, is 



I(P 2 



a 



16vr 2 



ft(21na- 1) + Vtt 2 ~ 1 



I n - vn 2 - i / 



+ 0(a 2 



(A.6) 
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where 



M± P 2 - mi c - Ml 

a = , it = 



2M±m Ac 



For P 2 — rrij^ ~ 0(p), which is the case for taking P 2 = mfi since ms c — tra c — 170 MeV, 
I (1) (P 2 ) starts from C(p 2 ), see Eq. (jBT6|) . The physical mass of the S c is above the A c 7r threshold, 
correspondingly ft > 1. 

Summing up the diagrams (a), (b) and (c), one gets the corrections to the mass of the S+ + at 

o( P 3 ) 

Am !oo P+ = - giM±^ + R e S (c) (mSe; M ± ,m Ae ), (A.7) 
where Re represents taking the real part. Similarly for the and S°, we have 



Am - T = -§^ + ^S( c )(m Ec ;M ,m Ac ), 



£2 

lot 



A loop a loop 



(Ai 



B Loop integrals 

Defining 



M 9 _ s - m 2 - M 2 

a = — , s = P z , ft = — 

m 2Mm 



A 



m 



(47T) 



' I(l n (4vr)+r'(l) + l) 



d-4 2 



the infrared singular part of the loop integral 

1 f d d k 1 

H ^ ~~i J (2vr) d (M 2 - & 2 ) [m 2 - (P - k) 2 } 

is (the expression for — 1 < ft < 1 is given in Ref. [44J ) 



Tf \ T( \ s-rr^ + M 2 
7(a) = 7(s) A, 



7( S ) 



1 



s 
a 



16vr 2 1 + 2aft + a 5 



[(ft + a)(21na-l) + F(s)] 



where 



Vft 2 - 1 lnf-v/^ 2 -!-^ 



In f Vft 2 - 1 - ft 



a 



(B.l) 



2\/l — ft 2 arccos 



v / l+2ar2+, 



7S / ' 



\/0 2 -l In Vfi 2 -l + fi + a -lnfi + a- VO 2 -! - 2ivr 



(B.2) 
< -1, 

-i< ft < i, (b.3) 
, ft > 1. 



In a more compact way, one can rewrite this by keeping the is, e — * + , explicitly, which is 
necessary to choose the correct Riemann sheet, 



F(s) = \/ft 2 - 1 In ( - yjQ 2 - I - Q, - a - ie) - \/ft 2 - 1 - ft - a + ie 



(B.4) 
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/W(s) is defined as 



1 f d d k W 



W t J (2ir) d (M 2 - k 2 ) [to 2 — (P — k) 2 ] ' { J 



One gets 



where in the IR method 

d d k 1 



A n = i I —— — ^-r — ^ — - — ^ = 2M 2 A H ^-lna , 

' (2n) d k 2 -M 2 V !6vr 2 / 

, d d k 1 , N 
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